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Abstract. 

This  paper  is  devoted  to  the  derivation  of  realistic  laws  relating  the  distribution  functions  of  reflected  and  incoming  particles 
at  the  wall.  A  polyatomic  scattering  kernel  is  derived  using  an  integral  operator  formalism.  Physically,  a  concept  of  partial 
accommodation  was  employed  in  the  kernel  derivation  in  agreement  with  different  accommodation  processes  of  the  three 
momentum  components  and  of  the  various  energy  modes. 


1.  INTRODUCTION 

Realistic  laws  relating,  at  the  wall,  the  distribution  functions  of  incoming  and  reflected  particles  are  needed  as  boundary 
condition  to  resolve  the  Boltzmann  equation  [1,  2,  3,  4].  Moreover,  in  the  moderately  rarefied  regime  these  laws  also 
allow  to  obtain  more  correct  velocity  slip  and  temperature  jumps  at  the  wall  and  to  broaden  the  validity  domain  of 
the  continuum  approaches.  In  this  topic,  a  polyatomic  kernel  is  derived  here  from  spectral  properties  of  its  associated 
integral  operator  using  a  well  known  scattering  kernel  formalism  [5,  6,  3]  The  molecules  are  considered  as  provided 
with  internal  structure  involving  rotational  and  vibrational  modes  and  are  possibly  in  a  strong  nonequilibrium  state 
concerning  the  vibrational  degree  of  freedom. 

The  eigenvalues  of  the  linear  integral  operator  so  considered  depend  on  five  parameters  appearing  in  the  kernel 
expression.  These  five  parameters  are  shown  to  be  equal  to  the  accommodation  coefficients  of  various  fluxes  at  the  wall, 
namely:  the  fluxes  of  the  three  components  of  momentum  and  the  fluxes  of  rotational  and  vibrational  energies  of  the 
molecules.  The  kernel  involves  a  part  obtained  in  another  papers  about  atom-like  molecules  [7]  and  so  also  allows  the 
modelling  of  anisotropic  behaviors  of  reflecting  surfaces.  Finally,  under  its  form  completely  developed  the  scattering 
kernel  appears  as  a  linear  combination  of  32  partial  kernels  describing  all  the  possible  associations  of,  respectively, 
diffuse  or  specular  processes  (according  to  three  directions)  and  elastic  or  inelastic  processes  (for  the  internal  modes). 
The  32  coefficients  of  this  combination  are  related  to  the  32  different  eigenvalues  of  the  integral  operator,  and  depend 
only  on  the  five  basic  accommodation  coefficients.  As  result,  in  addition  to  the  interplay  between  the  momentum 
components  included  in  the  kernel  of  unstructured  molecules,  the  polyatomic  model  also  allows  interplay  between 
the  different  energy  modes.  The  introduction  of  such  interplays  could  be  a  partial  reply  to  recent  criticisms  recently 
formulated  about  the  scattering  kernel  formalism  [8]. 


2.  POLYATOMIC  SCATTERING  KERNEL  DERIVATION. 

We  consider  the  problem  of  finding  the  scattering  kernel  B(V' ,  Eiri ,Eivi ,  giri ,  V,  Elr,  Elv.gjr)  governing  the  reflection 
of  polyatomic  molecules  at  the  wall.  V'  is  the  velocity  of  the  impinging  gas  particle  referred  to  the  wall,  V'  = 
( V'x ,  Vy ,  V'z )  €  { El!  =  El'x  x  El'y  x  Q!z  =  R  x  R  x  R}  and  V  the  velocity  of  the  reflected  one  referred  to  the  wall, 
V  =  (Vx.Vy,Vz)  €  {£2  =  Elx  x  £2V  x  O.  =  R+  x  R  x  R}.  These  velocities  reduce  to  the  peculiar  velocities  when  the 
slip  velocity  at  the  wall  is  neglected.  Vr  is  defined  as  Vr  =  (—Vx,Vy,Vz).  (x,y,z)  are  the  three  spatial  coordinates  with  x 
the  normal  axis  to  the  wall  oriented  from  the  wall  toward  the  gas.  Eiri  and  Elvr  are  respectively  the  rotational  energy  and 
the  vibrational  energy  of  an  incident  particle  at  the  wall.  Similarly  and  E,v  are  respectively  the  rotational  energy  and 
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the  vibrational  energy  of  a  reflected  particle  at  the  wall.  Then  the  subscripts  „  and  ,v  are  the  quantum  numbers  related 
to  the  internal  energy  of  a  particle,  hence  they  are  integers;  gjr  is  the  weight  of  the  rotational  degeneracy  and  will  be 
taken  here  equal  to  (2 ir+  1).  The  kernel  B,  which  is  the  density  of  probability  that  a  molecule  in  a  state  ( V' .  Eir/ .  Eivt ) 
hitting  the  wall  at  any  point  X  of  the  wall  is  reflected  at  the  same  point  in  a  state  (V,  Elri  Elv. ),  must  satisfy  the  following 
conditions  of  non-negativity,  normalization,  and  of  reciprocity  [3,  1,  2]: 


B(V  , Eiri , Eivf ,gir’,V, Ejr ,E/V, gir)  >0  , 


(1) 


V  f  B(V  ,  Eirf ,  Eiv!  ,  g,y ,  V,  Ejr ,  Ejy ,  gir)cly  —  1  , 

frz.J  n 


(2) 


Ill'll 


\V'\e  e-V e-V gir,B(V' ,Eir,,Eiv,,V,Eir,Eiv)  = 

__  l|v||2 

Vxe  e~Eire~Eiv girB(—V,Eir,Eiv,—V' ,Eitj ,Eivi)  ; 


(3) 


where  £,,  and  £,-,,  are  the  dimensionless  energies  :  £„  =  ,  £,,,  =  ,  with  k  the  Boltzmann  constant  and  Tw  the  wall 

temperature. 


2.1.  Analytical  derivation  from  integral  operator 


Let  us  write  the  transformation 


K(V,Eir, Eiv , gir, V' , Eirr , Eivt , giri )  =  [|  V'  |  MV\Eir,,Eiv,)]?[\  Vx  \  MV,Eir,Eiv)}^ 

E  ( V r  ■  EnJ ,  E[v/ ,  gjfj ,  V,Ejr,Ejv,gjr) 

where  fo(V.  Ejr,  Ejv)  is  the  local  equilibrium  distribution  function  at  the  temperature  Tw  of  the  wall,  defined  by: 

li  — Hv'll~  p .  „  £jr p  £,v 

f{)(V:Eir.Eiy)  =  ...  .  ^  8,r 


with 


Qr  —  y .  fi  ii 


(Cw0r)3 


QrQv 


a,  = 


2kTu. 


(4) 

(5) 

(6) 


Since  /o  is  a  known  function  the  problem  of  finding  B  is  equivalent  to  finding  K.  The  normalization  condition  and  the 
non-negativity  conditions  on  B  imply  obviously  the  same  conditions  on  K. 

Note  ,^r  the  set  of  the  rotational  energy  states  Elr,  and  ,':PV  the  set  of  the  vibrational  energy  states  Ely.  Consider  the 
five  elementary  Hilbert  space  of  states  L? (ElK)K=xyz,  Lr(^r),  and  L2(P?V)  of  square  summable  functions  with  their 
corresponding  usual  scalar  product 

<<Pku<Pk2>k  =  fnK-(pKi(V)(pK2(V)dv  ,  (pKi,(pK2  €  L2(klK)  ,  K  =  x,y,z 

<  (prl/(pr2  0 r  —  ^lir  *Prl  (^h'r)*Pr2(^h>)  ,  (p,  \  .  (px2  ^  E~  (,0 ) 

<  (jPy  |  .  (pv2  0y  =  tyvl  {Ejy)  (P\’2  {Ejy')  ,  G  L~  (-^ y) 

Consider  The  tensor  product  £  =  L2(£lx)  0  l2(kly)  ®L2(£lz)  0  L?  (00)  0'  L2  (.>0)  of  the  five  Hilbert  spaces  of  states, 
let  us  remark  that  this  tensor  product  of  Hilbert  space  £  is  dense  in  the  Hilbert  space  Jt?  =  L2  (£1)  ®  L2(/0r)  ®  Lr  (.>0) 
where  the  scalar  product  is  defined  by 

0  (p i  •  0—  y'.  /  t/h  fL,  Ejr,  Elv) (pi  ( V. ,  Ejr,  Eiv)dv  ,  (p i .  (p/  G  .  (7) 

jrjyJQ. 
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Instead  of  studying  the  problem  of  the  kernel  K,  we  study  the  linear  integral  associated  operator  A  defined  on  Jif  by 

MV)  =  T,  [  K(V,Elr,  Eiv ,  gir  ,V'.  Eir, ,  Eiv, ,  gir,  )y/(V'.Eir,,  Eiv, ,  glr,  )dv,  (8) 

,v  t,’  ■/£!' 


Assume  that  the  operator  A  has  a  purely  discret  spectrum,  and  assume  that  its  eigenfunctions  are  all  in  the  Hilbert 
space  £.  The  kernel  K  can  be  written  on  the  form: 


K=  L  %  Jy J,j.  Vi,  Vi,  (  l--ir )  Vn  ( i--n- )  (9) 

jx  ijy  ijz  dr  ijv  — 0 

Vjx  (V')y/jy  (Vy )  yjz  (V!)  y ir  (Eir, )  yjv  ( Eiv , ) 

where  the  functions  yjx (Vx)yj( Vv )Yj~(Vz) Wjr ( £)> )VjA Eiv )  are  the  eigenfunctions  of  A  with  their  corresponding 
eigenvalues  ^jx.jy.jz,jr,jv  ■  According  to  the  non-negativity  and  the  normalization  conditions,  the  eigenvalues  must 
satisfy  ^jx.jy,jz,jr,jv  G  [0,1]  for  all  jx,jy,jz,jnjv  G  N  .  Moreover,  one  can  see  that,  on  the  tensor  product  space  £, 
the  scalar  product  (7)  equals  the  scalar  product  defined  on  this  tensor  product  space  £  by  the  product  of  the  five 
elementary  scalar  products  (7).  So,  we  can  suppose  that  the  eigenvalues  have  the  form  ■'v  Av  and  that  the  set 

of  functions  yjx(Mx),  jx  €  N,  is  a  function  basis  of  the  %  corresponding  Hilbert  space,  X  =  x,y,z,r,v.  Therefore,  the 
expression  (9)  can  be  written  as  a  product  of  five  infinite  sums: 


K=  n  'L^hVjAMX)VhWx)  >  Mx  =Vx,Vy,Vz,Eir,Eiv  . 

Xe{x,y,z,r,v}  j= 0 


(10) 


Define  y0  =  VoxVovVozVorVov  by 


¥ox(Vx)  =  ^  |  Vx  | z  ,  y0v(Vy)  =  (CwVn)^e^  , 


C, 


- 1  p 

e  2  ,v 


¥0Z(VZ)  =  (C„V>0  ,  Vo AEir)  =  \j^e~'2e,r  >  VoAEn)  =  • 

Let  us  prove  that  t ]/q,  is  an  eigenfunction  of  A.  Mathematically,  the  normalisation  condition  can  be  also  written 

Y.  [  B{—  V, Eir,Ejv, gir,  —  VR,E,rf ,Ejvi ,giri)dvr  =  1  ,  (11) 

J  n' 


from  this  relation  (11),  the  reciprocity  relation  (3)  leads  to: 

L  [  I  V'  |  fo(V',Eir,,Eiv,)B(VAEir,,EiV,gir,,V,Eir,Eiv,gir)dv ,  =|  Vx  \  f0(-V,Eir,Eiv)  . 
jZTtj  Jn’ 


(12) 


Using  the  relation  (12),  the  calculation  of  A(y/o)  gives  A(yo)  =  yo  ■  Consequently  yo  =  VoxVoxVo-VorVov  is  an 
eigenfunction  of  the  operator  A  associated  to  the  eigenvalue  1 . 

Now,  following  the  five  microscopic  state  parameters,  let  us  introduce  five  parameters  related  to  the  eigenvalues  as 
follows:  Aoz  =  1,  and  for  j  A  0,  ^jx  —  (1  —  ax)  i°r  all  X  —  x,  y, z,  r.  v.  The  relation  (10)  becomes 


K=  FI  Wox{Mx)voAM'x)  +  ^-ax)1Lxi'h^Mx)VjAM'x)] 

X£  {x,y,z,r,v}  7=1 


(13) 


which  may  be  written 


K=  11  iaxVox  (Mx)vox  ( M'x )  +  (1 

Xe{x,y,z,r,v} 


ax)'LVh(Mx)Vjx(Mx)] 

7—0 


(14) 


Finally,  using  the  following  property 


L  V,AMx)V.iAMx)  =  S(Mx -M'x)  ■  (15) 

7=0 
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where  <5  is  the  dirac  function,  it  is  obtained 


K  =  {ax Yox(Vx)Y0x(V')  +  (1  -  ax)8(Vx  -  V')}{ay¥Qy(Vy)y0y(v;)  +  (1  -  0^)8  (Vy  -  Vv')} 
{azYoz(Vz) Voz(V')  +  (1  -  az)8{Vz  -  V')}{arY br{Vr)\i/0r{V')  +  (1  -  ar)8(Vr  -  V')} 

{«v  Y'ik (Vv)  Yov (Vy)  +  (1  -  av)8(Vv  -  Vv')}  • 


(16) 


Applying  inversely  the  transformation  (4),  the  operator  B  corresponding  to  the  kernel  K  above  (relationship  (16))  is 


?y  -Vt 

B(V' , Eir, , EiV , gir, ,  V, Eir,Eiv , gir)  =  {(1  -  ax)8(v;  +  Vx)  +  ax-^e  } 


(17) 


{(1  -  Oy)8(v;  -  Vy)  +  ocy-^^e  4  }{(1  -  ay)8(V'  - 
{(1  -  ar)8(Eir,  -  Eir)  +  arffe-£"-}{(  1  -  av)S{Eiv, 

In  the  further  calculations,  we  will  note  the  scattering  kernel  (17)  simply 


1 


Vz)  +  ocz 


-Eiv)  +  av-^-e  £,v}  . 


Cw\fn 

1 


B  =  PXPyPZPrPV 

where  Px.  Py.  Pz.  Pr,  Pv  corresponds  respectively  to  the  five  factors  of  the  expression  (17).  It  would  be  seen  that  these 
five  factors  satisfy 


/  ’  P,dVx  =  r  PydVy  =  /+°° PdVz  =  Y,Pr  =  1  •  (18) 

JO  J  OO  J  CO  -r  M 

On  the  other  hand  when  developing  the  expression  (17)  one  obtains  the  kernel  B  as  combination  of  32  elementary 
kernels  where  the  coefficients  are  functions  of  the  ax . 


2.2.  On  the  coefficient  ax 


In  this  section  we  prove  that  the  five  coefficients  ax  involved  in  the  scattering  kernel  equal  respectively  the 
accommodation  coefficients  of  the  various  fluxes  of  the  five  microscopic  state  parameters  ( Mx  =  Vx ,  Vy .  Vz .  Ejr ,  £IV ) . 
The  accommodation  coefficient  fi-x  of  a  physical  property  Mx  at  the  wall  is  defined  through  the  relation  [3?  ]: 


Px  ~ 


<&x  -o+ 

*x~*x 


(19) 


where  <t>x  is  the  incoming  flux  at  the  wall  of  the  property  Mx,  <P/  is  the  corresponding  reflected  flux,  and  <t>x  is  the 
reflected  flux  in  the  hypothetical  situation  of  perfect  accommodation  to  the  wall.  These  various  fluxes  are  written: 


cb 

X 


L  f  m\ V'x  Eir, , Eiv, , gir,)dv, 

ir'  ,iv'  ^ 


(20) 


= 

X 


E  [  m\Vx\Mxf+(V,Elr.Elv,glr)dv 


(21) 


where  /  and  f+  are  respectively  the  incident  and  the  reflected  distribution  functions  linked  by  the  relation 

|  Vx  I  fX  (y,Eir,gir,Eiv)  = 

E  [  I  vx  I  f~(y' ’Eir’,gir',Eivi)B(V' ,Eir,,girl,Eiv,,V,Eir,gir,Eiv)dV'  . 

WJa 


(22) 


974 


Accounting  for  (22),  the  reflected  flux  (expression  (21))may  be  rewritten: 

*$=  £  f  m\V'\  f-(V',Eir,,gir,,Eiv,)[Y,  f  MxPxPyPzPrPvdv\dv, 

ir'  ,iv'  ^  ^  ir,iv  ^  ^ 

and  the  reflected  flux  in  the  case  of  perfect  accommodation  is  written 

£  [  m\Vl\f-(V',E,y,gir,,ElV)[E  [  MxBedv]dv, 

ir.iv  ^ Q 


ct)e  — 

X 


where  Be  the  perfect  accommodation  scattering  kernel  is  defined  by 


BP  = 


2  gu 


QrQvCtn 


Vxe  cs-  e  £ire  e,v 


(23) 


(24) 


(25) 


2.2.1.  calculation  of  f3x 

The  normal  accommodation  coefficient  is  obtained  by  substituting  Mx  ='[  Vx  |  in  the  definition  (19).  In  this  case, 
accounting  for  the  property  (18)  and  the  expression  of  the  partial  operator  Px  it  is  obtained 

£  /  |  Vx  |  PxPyPzPrPvdv  =  -(1  -  ax)V'  +  axCw/7Z 
777.,  J  Cl  2 


then  the  expression  (23)  yields 


<£>;  —  4>  =  ax 


,£  [  m\V;\f-(V\Eir,,Eiv,,gi,J)(V'  +  C^^)dv,  . 
if%SJof  2 


The  calculation  of  (t>r  leads  easily  to 


/  Cw  yTr 


[  m\V'\f-(V',Eir,,Eiv,,gir,)(V'+^X^)dv,  . 


Consequently  we  obtained  from  the  relation  (19)  j 8X  =  ax 


2.2.2.  calculation  of  /3V  and  pz 

The  tangential  accommodation  coefficient,  /3V  is  obtained  by  substituting  Mx  =  Vy  in  the  definition  relation  (19).  In 
this  case  it  is  easily  seen  that  <Pev  =  0.  Then  accounting  for  the  property  (18)  the  calculation  of  expression  (23)  leads  to 
Py  =  (Xy.  Similarly  it  is  found  /3,  =  a-. 


2.2.3.  calculation  of  p,  and  j8v 


Now  substitute  Mx  =  girEjr  in  the  relation  (19).  Accounting  for  the  property  (18)  it  is  obtained 


£  [  girEirPxPyPzPrPvdv  =  Y^girEirPr 

ir.iv'^  ir 

where  we  have  noted 


(1  ar)gjriEitj  +  a, 


q; 

Qr 


q;  =  Y^slEi, 


ir 
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Then  the  expression  of  <£>+  leads  to 


=  a, 


L  [  m  |  V'x  \f~  {V ,  Eir, ,  Ew ,  girf )  (g/r'Ei, 

;7 J  O! 


-I'"'" 


Using  the  expression  (25)  of  Be ,  we  obtain 


^  /  girEirBedy 

irivJQ. 


Q £ 

Qr 


and  then 


L  /  /n|TA,|/-(V,,£,y,JE,v,g,y)(g,yZi,-; 


-I’*' 


consequently,  /3,  =  a,-.  In  the  same  way,  substituting  =  £)v,  it  is  found  /3V  =  av  .  In  conclusion,  the  five  parameters 
involved  in  the  scattering  kernel  (17)  are  the  accommodation  coefficient  corresponding  to  the  five  state  parameters, 
namely  the  three  momentum  components  and  the  two  internal  energy  degrees. 


3.  CONCLUDING  REMARKS 

Under  its  factorized  form  (17),  the  proposed  scattering  kernel  is  easy  to  use  in  the  analytical  calculations  or  to 
implement  in  numerical  modelling,  notably  because  of  its  simple  factorized  form  where  the  incoming  molecule 
parameters  appear  only  through  Dirac  functions.  In  order  to  show  its  physical  meaning,  the  expression  (17)  may  be 
developed.  Under  its  developed  form,  the  scattering  kernel  appears  as  a  linear  combination  of  32  elementary  scattering 
kernels.  All  these  elementary  kernels  correspond  to  various  situations  of  accommodation  at  the  wall  [7],  The  linear 
combination  coefficients  which  represent  the  weight  of  the  various  types  of  accommodation  in  the  reflection  process 
are  combinations  of  the  factors  ax  and  (1  —  otx).  In  each  elementary  kernel  each  molecules  states  accommodates 
independently  from  the  others.  So  the  new  kernel  allows  to  take  into  account  the  interplay  between  the  molecule 
freedom  degrees  when  interacting  at  the  wall  [8],  Finally,  the  proposed  scattering  kernel  appears  as  a  complete 
construction  based  on  the  five  accommodation  coefficients  corresponding  to  the  five  microscopic  parameters  defining 
the  state  of  the  molecules. 
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